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Abstract 

The exact relativistic form for the beta decay endpoint spectrum is derived and presented in 
a simple factorized form. We show that our exact formula can be well approximated to yield 
the endpoint form used in the fit method of the KATRIN collaboration. We also discuss the 
three neutrino case and how information from neutrino oscillation experiments may be useful in 
analyzing future beta decay endpoint experiments. 
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I. INTRODUCTION 

The di=cove.y ot neutrino oscillation. Q fl fl 3 prober the neutrino squared n,ass 

differences and mixing angles p, but leaves open the issue of what is the absolute scale of 
neutrino mass. The latter has important cosmological implications in the cosmic microwave 
background and large scale structure in the Universe, as already indicated by the sensitivities 



reached, for example, by the recent WMAP-3 j?], the 2dF Galaxy Redshift Survey 
Sloan Digital Sky Survey results [91]. One expects even better sensitivities in the next 



and 



generation of cosmological observations 10|, lUl]. Interesting as these may be, there are 



essentially only two ways to get insight into the absolute scale of neutrino mass in the 
laboratory: searches for neutrinoless double beta decay [12|] and investigations of the beta 
spectra near their endpoints [isl . Q, Q, Q, 0, [isl . For the latter direct search for the 
neutrino mass a very low beta endpoint is crucial: tritium was used in the most sensitive 
spectrometer experiments [l4 , 
19|. 



15l | and rhenium in the up-coming cryobolometer experiments 



Currently a next generation tritium beta-decay experiment is being prepared, scaling up 
the size and precision of previous experiments by an order of magnitude, and increasing 
the intensity of the tritium beta source: the KArlruhe TRItium Neutrino experiment KA- 
TRIN 13, IgI, 0, Such an improved sensitivity experiment will probe neutrino masses 
ten times smaller than the current limits and therefore play a crucial role in probing for 
direct effects of neutrino masses. 

Prompted by the prospects that high sensitivities can be achieved in the next generation 
of high precision neutrino mass searches from tritium beta decay experiments la, llTI] we 
reexamine the accuracy of the kinematical formulae used in the determination of neutrino 
masses from the shape of the endpoint spectrum. We also discuss the interplay of neutrino 
oscillation data and the expectations for the beta decay endpoint counting rates for the 
different types of neutrino mass spectra. 

II. RELATIVISTIC BETA DECAY KINEMATICS 

In what follows we label the relativistic momenta and energies involved in tritium beta 
decay according to 



(1) 



2 



J E E E' ^ ' ^^'"""''"'^ -P - Pe-Pu), (3) 



The masses of ^He'^, e~ and Ue are denoted by M', rUe and respectively. In order to see 
the convenience of an exact relativistic description we mention, as recently noted in Ref. [2^ 
that the well known relativistic formula for the maximum electron energy 

= ^[M' + ml - (m. + M'f] , (2) 

gives a value about 3.4 eV lower than the approximation M — M' — ruy often used. This 
suggests the desirability of carrying out the full phase space integration using relativistic 
kinematics. 

Start from the standard formula for the decay width at rest, 
^ If d^ped^pud^p\ 

where denotes the spin-summed, Lorentz invariant "squared" amplitude. To explore 

the constraint of Lorentz invariance one might a priori consider expanding in terms of 

invariants constructed from the four-momenta. For example, up to two powers of momenta, 
the most general form is, 

\M\^ = A- Bpe ■ Pu - Cp ■ Pinitial + (4) 

where A, B and C are constants. Now it is easy to perform some initial integrations. As 
usual J d^p' is first done with the momentum delta-function. Then the angle between pe and 
Pi, is eliminated using the energy delta function. Three more angular integrals are trivial. 
As the result one may replace in the '^H rest frame 

\M\^ -^A + B{E,E, - Pe • p.) + CM{M -E,- E,), (5) 

where 

Pe ■ p. = ^ [M^ - M'2 + ml + ml- 2ME, + 2E,{E, - M)] . (6) 



Eq. can now be inserted in the resulting usual formula |2l| 

1 



dE^dEe\M\'. (7) 



Next we find dT/dEf. by integrating over dE^y for each Ef.. The limits of integration E^^^{Ee 



and E^°'^{Ef,) can be read from 2l|. The most tedious part of the present calculation is 
finding the factorizations: 

2m, M'' 



rpmax rpmin ■^-'"/^e / rpmax rp \l/2 



M 



2M 
wherein: 

(mia)^ = - 2MEe + ml (10) 



The importance of the factorization is that it makes the behavior at the endpoint E^, = E\ 
transparent. Then we have the exact relativistic resuh, 



max 
e 



'y{y + + CM(M - Ee) + 



where ?/ = E^'^^ — E^. 

As it stands, this formula is based only on the kinematical assumption in Eq. (j4j). It 
obviously vanishes at the endpoint ?/ = as ^Jy. Note that all other terms are finite at 



y = 0. The overall factor \/y{y + IrriyM' jM^ gives the behavior of extremely close to 
y = for any choice of A, B and C, but departs from away from the endpoint. 



Dynamics is traditionally put into the picture [22] by examining the spin sum for a 
4-fermion interaction wherein the nuclear matrix element is assumed constant. This is 
presented as a non- Lorentz invariant term, 

\M\'' = BEeE,. (12) 

We will see that this is excellently approximated in our fully relativistic model by, 

A = C = 0, B^O. (13) 

A more accurate treatment of the underlying interaction might give rise to small admixtures 
of non-zero A and C as well as other unwritten coefficients in Eq. (jlj) above. 

The form for the spectrum shape near the endpoint that results from putting A = C = 
in Eq. is 



dEe (27r)34(mi2)4' e^^^y^' ^ 



{MEe-mi)Jy[y + ^— \ y + ^{M' + m,) . (14) 



M 



Note that if we had employed the non-relativistic form given in Eq. f|T2|) the net result would 
be a replacement of an overall factor in Eq. f|T^ according to. 



{ME,-ml)-^{ME,-E'^). (15) 



The difference of these two factors yields the contribution of the pe • Pu term. It is really 
negligible near the endpoint region since it is proportional to pi and is suppressed like 
pl/{MEe) compared to unity. We have checked that the result of our calculation with just 
the EeEi, term agrees with the calculation of Ref . , though their result looks much more 
complicated, as they did not present it in the simpler factorized form given here. 

Note that only the two rightmost factors vary appreciably near the endpoint of Eq. (|Til) . 
If we further approximate M'/M 1 and ^^'^J^" — > 1 the endpoint shape is well described 
by 

dT 



dE. 



a {y + my)^Jy{y + 2my). (16) 



Now we compare with the formula used in the experimental analysis 15 1 

dV 



a 



{Eo -V,- E)^{Eo -V,- Ef - ml (17) 



dE, 

This agrees with the above approximation in Eq. ( flGil if one identifies 

{Eo-Vi-E)=y + m,. (18) 

Note that E is the non-relativistic energy given hj E = E^ — me- Furthermore, Eo — Vi is 
identified with our (M - M' - - 5^™"^). dE'^"'-' is defined by, 

^max _ - M' - - dE'^"'', (19) 

and was shown in 20| to be independent of m^ to a good approximation. Thus we see that 
the exact relativistic endpoint structure obtained here may be well approximated by the 
form used in the experimental analysis. 

Often, authors express results in terms of a variable, x, which from our discussion may 
be seen to be the same as, 

x = -y-m, = E,- E,"^"^ - m,. (20) 

In Fig.lll dr/dEe as computed from the exact formula, Eq. ( fT4l) is compared with its approx- 
imate analog as a function of x. As can be seen, the differences between the approximate 
and exact formulae are tiny. 

It may be worthwhile to remark that the exact relativistic kinematical expression in 
Eq. (|T4l) is no more complicated than the approximation one ordinarily uses. 
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Figure 1: Comparison of approximate and exact formula 



III. NUMERICAL SIMULATIONS 



We have done Monte Carlo simulations by creating random data sets following the exact 
relativistic kinematical expression in Eq. f|T^ for neutrino masses of = eV and mj, = 
1 eV and fitting them with the standard formula Eq. f[T7|) . The simulations were performed 
for a KATRIN-like experiment {l8| considering: 



The ro- vibrational states of the T2 {T^He)^ decay |2J] 



• A signal rate from a KATRIN-like molecular gaseous windowless tritium source with a 
column density of 5 ■ 10^^ molecules/ cm^ over an active area of 53 cm^ and an accepted 
sohd angle of AVL/At: = 0.18 

• An expected background rate of 0.01 s~^. 

• A response function of a KATRIN-like experiment considering the energy losses within 
the tritium source and the main spectrometer transmission function with a total width 
of 0.93 eV. 

• 3 years of total data taking covering an energy range of the 25 eV below and 5 eV 
above the tritium endpoint following an optimized measurement point distribution 

H. 

Fig. [2] shows the results for the observable obtained from the fitting of 10000 sets 
of Monte Carlo data randomized according to the exact relativistic formula Eq. ( fT4l) and 
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Figure 2: Results on the observable mf, from the fitting of 10000 sets of Monte Carlo data ran- 
domized according to the exact relativistic formula Eq. ()14|) and to a fitting routine based on the 
standard formula Eq. (jl7p . for neutrino masses of eV (left) and 1 eV (right). The rms values 
of the Gaussian-like distributions correspond to the expected statistical uncertainty Am^ gf^f for a 
KATRIN-like experiment. 

to a fitting routine using the standard formula Eq. (|T7j) . assuming neutrino masses of eV 
(left) and 1 eV (right). The rms values of the Gaussian-like distributions correspond to 
the expected statistical uncertainty Am^^^^^ for a KATRIN-like experiment. Clearly the 
mean value of the fit results for the neutrino mass squared does not show any significant 
deviation from the starting assumption of = eV or = 1 eV, respectively. This 
establishes that the exact relativistic formula Eq. ( JT^ can be well approximated by the 
standard equation f|T7|) for the precision needed for the next generation tritium experiment 
KATRIN. This is probably due to the fact that KATRIN is investigating the last 25 eV 
below of the beta spectrum below its endpoint only, where the recoil corrections are nearly 
independent on the electron energy. 
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IV. THREE NEUTRINO CASE 



Of course, the most interesting application is to the case of three neutrinos with different 
masses, mi, m2 and ma. Then there will be a different endpoint energy, E^^"'^ corresponding 
to each one. The effective endpoint factor in the good approximation of Eq. (fT6l) is the 
weighted sum, 

3 

F,jf{E,) = J2 iKuWvi + m,)[y,{yi + 2m,)f^e{y,), (21) 



i=l 



where yi{E^) = E^"-^ — E^ and the Ku are the elements of the 3x3 lepton mixing matrix 25|, 



261 ]. We note that the further good approximation that the quantity SE^""^ is independent 



of the neutrino mass, gives the useful relation 

yi - yj = rrij - rui. (22) 

Now let the unindexed quantity y stand for the yi with the smallest of the neutrino masses. 
Using Eq. fl2^ allows us to write the explicit formula for the case (denoted "normal hierar- 
chy") where mi is the lightest of the three neutrino masses as: 

F^niy) = \Kn\\y + m^)[y{y + 2rm)]'/' (23) 
+ \Ki2\'^{y + mi)[(?/ + mi - m2)(y + mi + m2)]^^'^6{y + rrii - 7712) 
+ \Ki3\'^{y + mi)[{y + mi - m3)(?/ + mi + ms)Y/'^9{y + mi - mg). 

In the other case of interest (denoted "inverse hierarchy" ) we have: 

FjHiy) = \Kn\\y + msMy + 2ms)Y/^ (24) 
+ \Kii\^{y + m3)[{y + ms - mi)(y + ma + mi)Y^^9{y + ms - mi) 
+ \Ki2\^{y + ms)[{y + mg - m2)iy + mg + m2)Y^^9{y + mg - ms). 

where mg is the lightest of the three neutrino masses. From these equations we may easily 
find the counting rate in the energy range from the appropriate endpoint up to ymax as 
proportional to the integral 

/Vmax 
dyFNHijj), (25) 


or, for the "inverse hierarchy" proportional to, 

/Vmax 
dyFiuiy). (26) 




We note that, as stressed in ref. [20], information on neutrino masses and mixings ob- 
tained from neutrino oscillation experiments is actually sufficient in principle to predict 
n{ymax) as a function of a single parameter (up to a twofold ambiguity). Thus, in princi- 
ple, suitably comparing the predicted values of n{ymax) with results from a future endpoint 
experiment may end up determining three neutrino masses. 

To see how this might work out we make an initial estimate using the best fit values 
of neutrino squared mass differences, 

(27) 



A = 


2 2 
1712 ~ '' 


= 7.9 X 10-^e\/2^ 


B = 


1 2 2 

mg — m2 


= 2.6 X 


coefficients. 










= 0.67, 






= 0.29, 




I^13P 


= 0.04. 



(28) 



Currently l-ft'is] is consistent with zero and is only bounded. For definiteness we have 
taken a value close to the present upper bound. However, we have checked that the effect 
of putting it to zero is very small. Now, from the two known differences in Eq. ( 1271) we 
can for each choice of (considered as our free parameter) find the masses mi and m2, 
subject to the ambiguity as to whether is the largest (NH) or the smallest (IH) of the 
three neutrino masses. Of course we hope that future long baseline neutrino oscillation 
experiments 27, 28, 0, [s^ might eventually determine whether nature prefers the NH or 



the IH scenario. 

Fig. [3] shows typical solutions for the mass set (1711,1712) in terms of the free parameter 
7713. Very large values of would fall within the sensitivity of upcoming cosmological 



tests 
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In the right panel of Fig. [3] we display the predicted values of n{ymax) for 
each possible mass scenario and the choices of (1,10,20) eV for Umax- These quantities are 
proportional to the electron counting rate in the energy interval from the endpoint (for 
each mass scenario) to Umax eV below the endpoint. The different values of Umax reflect, 
of course, different experimental sensitivities. The main point is that, for sufficiently large 
m3 values, the counting rate is seen to distinguish the different possible neutrino mass sets 
from each other. We hope that the present method of relating observed neutrino oscillation 
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Figure 3: The left panel shows typical solutions for (mi, 771,2) as a function of 7713 for the NH case 
(solid curves) and the IH case (dashed curves); the middle dot-dashed is given for orientation. The 
right panel give the predictions for the quantities, n{ymax)i proportional to the event counting rate 
which includes emitted electrons within, respectively 1 eV, 10 eV and 20 eV from the appropriate 
endpoint. 

parameters to predictions for the beta decay endpoint counting rates may play a useful role 
in the forthcoming experiments. 

V. SUMMARY AND DISCUSSION 

We have derived the exact relativistic form for the beta decay endpoint spectrum and 
presented it in a very simple and useful factorized form. We showed that our exact formula 
can be well approximated to yield the endpoint form used in the fit method of the KATRIN 
collaboration. This was explicitly established through a detailed numerical simulation. We 
have also discussed the three neutrino case and shown how information from neutrino oscil- 
lation experiments may be useful in analyzing future beta decay endpoint experiments. 
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